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Abstract 

A recurrence relation of Riccati-type differential equations known in super- 
symmetric quantum mechanics is investigated to find exactly solvable poten¬ 
tials. Taking some simple ansatze , we find new classes of solvable potentials 
as well as reproducing the known shape-invariant ones. 
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It is generally difficult to solve exactly the eigenvalue problem of a (time-independent) 
hamiltonian in quantum mechanics. Among the various methods developed for this pur¬ 
pose, a remarkably simple but powerful one is a ladder operator technique, whose typical 
example is a simple harmonic oscillator. If a hamiltonian has a discrete eigenvalue spectrum 
bounded from below, the energy eigenstates should be labeled by integers and formal rasing 
and lowering operators can be written in this basis. However, it does not provide a way to 
find the ladder operators explicitly for a given hamiltonian. A practical method of getting 
ladder operators has been studied by several authors( ||- ||) using ideas of supersymmetric 
quantum mechanics( |l|], [||) and a concept of shape-invariant potentials [[|. This approach 
has (re)produced many exactly solvable potentials. In this paper, we investigate the recur¬ 
rence relation of Riccati-type differential equations arising in the supersymmetric quantum 
mechanics directly taking some ansatze. In addition to reproducing known shape-invariant 
potentials concisely, new classes of solvable potentials are found. 

Let us begin by summarizing some results of supersymmetric quantum mechanics which 
are used for the present work(for a review, see [2J for example). Consider a family of one 
dimensional hamiltonians 




( 1 ) 


where p = —, which are connected by 


Hk+i — B k B\ — E k ( 1 ) — B\ +l B k+ i , 



( 2 ) 


Here B k and B k are defined using superpotentials a k (x) 



( 3 ) 


where a k (x) is a real function such that the ground state wavefunction which is proportional 
to exp(— f a k (x)dx ) should be normalizable. The nth eigenvalue of H k is denoted by E k {n ) 
n — 0,1, 2, ■ • •, and a relation between eigenvalues is given by 


E k +i(ji) — E k (n + 1) — E k ( 1) — Eq(ji + k + 1) — E${k + 1) . 


( 4 ) 
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Note that all the ground state energies are defined to be zero. The E k ( 1) is assumed to be 
positive and is also denoted by s k , standing for a level spacing between the kth and (/c + l)th 
eigenstate of the starting hamiltonian H 0 . It is worth noting that E 0 (n + 1) = YJk =o s k , n = 
0,1,2,- • ■ and ifo(0) = 0 . Denoting the nth eigenfunction of H k as r ijj k ^ n )(x), important 
relations between the eigenfunctions follow from (j2j) as 

i>k+ l(n)O) = {E k {n + 1 )}^ 1/2 Bki>k(n+l){x ) , 

B ] k ^ k +l(n)( x ) = {E k {n + l)} 1/2 ^fc(n+l)(^) • (5) 

Note that B k i/) k (o) (x) — 0 . A repeated use of this relation gives 

-Bo-Si • • • B l i>k+i{n){x) = {E 0 (n + k + l)E 1 (n +k) ■ ■ ■ E k (n+1)} 1/2 ^ n+k+1 )(x) . (6) 

Putting n = 0, this equation can be used to get the excited-state wavefunctions of H 0 from 
the knowledge of the ground states of systems related as superpartners. 

Now, putting B k and B k in ([]) into ([|), a recurrence relation of Riccati-type differential 
equations is obtained as 

2I4 + i (x) = a 2 k+1 (x) - a' k+1 (x) = a 2 k (x) + a' k (x) - 2s k , (7) 

where the prime means a derivative with respect to x . To search for the possible family of 
a k (x) and s k satisfying the difference differential equation (0), it is convenient to put a k (x) 
and a k +i(x) as 

«fcO) = 9 » 1 / A 2s k f k (x) +f k (x)} , afc+i(a;) = {2 s k f k (x) - f£(x)} , (8) 

z J k \ x > z J k \ x ) 

which solve (0) identically. It holds for non-zero s k only, and in general, there can be added a 

free parameter a k , say, inside the curly brackets in (||). Note that f k (x) is proportional to the 

first excited-state wavefunction divided by that of the ground state of the fcth hamiltonian. 

The consistency in the two equations of (||) requires 

{fk( x )f k +i( x )Y = 2 {s k f k ( x )f k+ i( x ) ~ s k+if k ( x )fk+i(x)} ■ (9) 
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We further rewrite equations (||), © by putting f k {x) as following : 

fk(x ) = g k (x) exp(J h k (x)dx) . ( 10 ) 

Note that there is an ambiguity in the choice of h k (x), or equivalently in g k (x) : Adopting 
g k (x ) = g k (x) exp (f /3 k (x)dx ) , and h k (x ) = h k (x ) — /3 k (x ), instead of tilded functions, f k (x ) 
remains the same for arbitrary (3 k {x). We remove this ambiguity by choosing j3 k {x) such 
that ho(x) = 0 , and 

fk(x ) = g k {x) exp (J h k {x)dx) , f' k (x ) = fo(x)exp(J h k {x)dx) . ( 11 ) 

The consistency in this procedure requires for all k = 0,1, 2, • • • , 

/o(®) = + 9k(x) h k (x) , (12) 

then the difference differential equation (0) is satisfied. Now, equation (^) becomes 

fo( x ) + ^ fo(x) ( h k {x) + h k+1 (x)) = s k g k (x ) - s k+1 g k+1 (x) . (13) 

Solving this equation for g k (x) and inserting it into (p~2|) yields 

kfo(x) + p k (x)fo(x) + (q k (x) + s k - s 0 ) f' Q (x) - s 0 h k (x)f 0 (x ) = 0, (14) 

where 


Pk{x) = khk(x) + 9 k (x) , q k (x) = 0' k [x) + h t (x)6 t (x) , 


(15) 


and 


fc-i 


0 k (x) = - ( h j( x ) + kj+i( x ) ) , 0 o (x) = 0 . 


(16) 


3=0 


Since (|14[) is identically satisfied for k = 0, it is equivalent to the difference of (k + l)th and 
kth expressions in (|T4j). Thus, with the difference denoted by A y k = y k+1 — y k , we finally 
obtain 


fo(x) + Ap k (x)fo(x) + (Aq k (x) + As k ) fo(x) - s 0 Ah k (x)f 0 (x) = 0. (17) 
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It is equivalent to the difference differential equation (0), and we will try to find fo(x) and 
Sfc consistently by making simple ansatze on h k {x) below. The superpotential in (|j) is now 
given as 

M x ) = y 77 ^y{s°/o(^) - (k-^)fo(x)} + i/i fc (x) - 6 k {x) . (18) 

Note that an overall constant factor of fo(x) does not affect the superpotential. 

Let us solve ([17]). The simplest ansatz on h k (x) may be putting h k (x ) = 0 for all 
k — 0,1, 2, - - • . In this case, ([[7]) becomes 

/o" 0*0 + ^Skfo(x) = 0. (19) 

This is solved consistently by putting 


As k = s, 


( 20 ) 


where s is a constant parameter. It is easy to find fo(x) in (|1^), and the superpotential is 


given in ([18]) with h k (x) — 0, 6 k (x) = 0 in this ansatz , that is, 


M x ) = -Frrr{ s ofo( x ) ~ (k - ]-) ft(x) } . 
Jo\ x > z 


( 21 ) 


The eigenvalues are given by the common formula 


E k (n) = s 0 n + - sn(n — 1 + 2k) , 


( 22 ) 


referring to (20) and (|l|). In the following, we list the resulting superpotentials, while the 
potentials can be written down explicitly by using 2V k (x) = a\(x) — cx' k (x ) in (|^). With 
constant parameters 6, d , and putting 

2( Sq T ks ) — s 


V = a k = 

(1) s>0 ; fo(x) = sinrj(x + d) 

—a k cotrj(x + d) + 


2 r] 


a k (x) = 


tan % (x + d) 


sinrj(x+d) 

a k~ b r-rX R 


Eckart potential 


(23) 


cot | (x + d) Poschl-Teller potential I 
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(2) s<0 ; fo(x) = c\ sink r]{x + d) + c 2 cosh rj(x + d) with constant parameters c \, c 2 . 
Thus, the superpotential in (P is 

b 


a k ( Ci + c 2 tank r](x + d )) 
OLk\x) = -:-:-t—t—— x-h 


c 2 + ci tanhrj(x + d) ci sinhr](x + d) + c^coshpix + rf) 


Some well-known potentials in the literature come out : 


(24) 


a. /o(x) = sinhr}(x + d) (ci = 1, c 2 = 0 ) 


a fc (x) = 


Rosen-Morse potential 


a k cothri(x + d ) + 

tank \ [x + d) + 4^ co ^ 77 (x + d) Poschl-Tcllcr potential II 


(25) 


b. Morse potential ; /q(x) = coshpix + d) ( c\ = 0 , c 2 = 1) 

b 


a k {x) = a k tanhrj(x + d) + 


coshr}(x + d) 


(26) 


c. Morse potential ; /g(x) = expp(x + d) , (ci = c 2 = |) [ ci = —c 2 = — ~ also gives a 
similar potential.] 


cx fc (a;) = a k + bexp(-r/(x + d)) . 


(27) 


(3) s = 0 ; two cases arise as following, 
a. Harmonic oscillator potential ; /q(x) = 1 

a k (x) = s 0 (x + d) . (28) 


b. Harmonic plus inverse-square potential ; /q(x) = x + d 

a k( x ) = ^-(x + d) -, X k = k - i + b . (29) 

2 x + d 2 

It should be noted that k has a maximum allowed value when s < 0 since s k = So + ks is 
assumed to be positive. It is remarkable that a large part of the known exactly solvable 
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potentials given in refs. 0, 0 are reproduced here within a single ansatz. Further details 
of the potentials and the wavefunctions can be found in refs. §, 0, §, 0 where a notion 
of shape-invariance of potentials is utilized in the framework of supersymmetric quantum 
mechanics. 

Let us solve (jT7|) with an ansatz which is slightly more general than assuming h^{x) = 0 
for all k , 

! 0 for k = even 

(30) 

h(x ) for k = odd . 

Using for k = even, odd integers respectively, 

A h k {x) = ±h(x) , A p k (x) = h(x) {1 ± (k + } , A q k (x) = {1 ± (2k + 1) } , 


the equation (|T7|) reduces to 


fo( x ) + ^h( x )fo( x ) + { + 

/"w - \hx)K(x) + _ 



+ As 0 } fo(x) - s 0 h(x)f 0 (x) = 0 , 
+ As i } fo( x ) + s 0 h(x)f 0 (x ) = 0 , 


(31) 


for k — 0,1 and 

/ » W + { ¥ + 2SR> / »W = 0 ' (32) 

which results from the consistency of (|l^) for all k = 2,3,4,5, - • •, together with A Sk 
determined with a constant w as 


A s 2 k = As 0 + kw , As 2 fc+i = Asi — kw 


(33) 


Combining equations (^l|) and (|32|), one finally arrives at 


h '( x ) + \h 2 (x) (- + 1) + Aa = 0 , 2af' Q {x) = h(x)f 0 (x) , 
z a 


w = 2A a , (34) 


where A, a are 


A — Asq Asi , 


a = 


Aso — Asi — w 
4 Sn 


(35) 
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We consider a^O, since otherwise h(x) = 0 in (|34|) , which has been already covered. The 
energy level spacings are now given as 


S 2 k — Sq + k\ , S 2 k +1 — (1 + 2cr) { s 0 + (k + —) A } . 


(36) 


Here, cr > —1/2 for positive level spacings. From this and (4j), the energy levels are given 


as 


E k {n) = (1 + a) n { s 0 + j (n - 1 + 2k) } 

l-(-l) r 


A 


A 71 


(-1)M ~ 2 + j(2*-l)) - '-f (-1)"}. 


(37) 


This formula reduces to (|2^) as cr —^ 0 , by noting in ([56]), A Sk = A/2 in this case. The 
superpotential in ( |l8|) is given finally using (|34p 

2 o , A Xk N h{x) 


to(x) = y) - (-It 


(38) 


h(x ) v ~ u 4 ' 2 y v ' 4 

Now, we only need to find h(x ) satisfying ([54]). The resulting superpotentials are listed 
in the following with a constant d , and putting 

a(-l) fc 


V = 


|A| ( 1 + cr) 


i / t w A A k . 

h = (1 + <r) (so - 4 + y). to = 2(1+ a) 


(1) A > 0 ; h(x) = cot tj(x + d) 


tan rj (x + d) 

a k {x) = b k - jd k r]cotr]{x + d) . (39) 


V 


(2) A<0 ; h(x) = |f -cothrj(x + d) 


. . , tank rj (x + d) . 

a k {x) = b k - — /Ufc rj coin q [x + d) . 


71 


(40) 


(3) A = 0 -h(x) = 


cr fc (a;) = s 0 (1 + a) (x + d) 


x + d 


(41) 
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These potentials are similar but different from the ones obtained by our first ansatz , that 
is, Poschl-Teller potentials I, II given in (^3|), (|25|) respectively, and (|29|). This fact is clear 
in view of their different energy-level structures and the wavefunctions. For example, the 
potentials corresponding to the superpotentials in (^), ([14]) are given respectively (putting 
d = 0 for simplicity) 


V k (x) 


-^-x‘+ 


Ag-A fc 

2x 2 


so 

2 


1 N 

(Afc + -j), 


(42) 


v k(x) — ~ s l(l + a) 2 x 2 + 


|(l + «T+(-l) t «T). 


(43) 


The energy level spacings are s k = Sq for all k, for the potentials (fEj), while s k = So , So (1 + 
2cr) for k = even, odd integers respectively for (^) . The difference between the two systems 
can be seen more clearly in their wavefunctions. The unnormalized wavefunctions for Vq{x) 
given in ([42D are 


/ ( \ Ao / I 2 \ 7-Ao — v 

Wo(k){x) = x°exp{--s 0 x )L k ( ) , 


while for Vq(x) given in (43) 


(44) 


^ 0(2 k)(x) = x^°exp(~-s 0 {l + a)x 2 ) L k 2 (s 0 (l + a)x 2 ) , 


i> 0 ( 2 k+i)(x)=x exp(--s 0 {l +a)x )Ll (s 0 (l + cr)ar) 


(45) 


where L“(^) is a generalized Laguerre polynomial. Note that as a —> 0, the potentials in ([43|), 
the wavefunctions in (^) and the energy levels all reduce to those of the simple harmonic 
oscillator. In contrast, if the parameter is chosen such that the potentials in (|42|) become 
quadratic, only ‘half’ of the energy levels and wavefunctions of a corresponding harmonic 
oscillator are obtained as can be seen in (|44|). We remark that the wavefunctions in ([45]) 
are normalizable as long as a > —1/2, although their first derivatives are not continuous at 
x = 0 . The point is that another class of potentials like in ( |43| ) is possible as well as the well- 
known ones like in ([E]) if we seek supersymmetry-related potentials possessing normalizable 
wavefunctions. 
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How can the (radial) Coulomb potential consisting of 1/x and 1/x 2 terms be reproduced 
in our approach? Since we need h k (x) as inputs for solving (|17|), we assume h k {x) to be that 
of Coulomb potential 


hk(x) = g ( 


) 
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(46) 


k -\~ CL k \ CL Qj(cL T 1) 

where g, a are constant parameters. The derivation of this will be described later together 
with other general idea of getting superpotentials. Then ([[7]) becomes for k = 0, 

9 + b }fo( x ) ~ soh f 0 (x) = 0, 

(47) 


fo( x ) d- TT fo ( x ) + { ( ^ + 9 ^, y ) hi + 1 x 2 

2 a(a +1) a z (a + l) z 


where hi = — 2g/{ a(a + l)(a + 2) } , and 

(2o+l)/ 0 » - {2A + 2g(2 " +1) }&(*) + 2 soMx) = 0, (48) 

a(a + 1) 

which results from the consistency of (p~7|) for k — 1, 2, 3, • • • . Here, s k is determined as 

.2 1 i 


Sfc 
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) + \h k + bk + c , 


(49) 


2 (k + a) 2 (k + 1 + a) 2 
with another constants A ,b,c. Combining ([47l) and (48), interesting cases arise only when 


1 1 

A = 0, c = (o+ —) 6 , ci=t -- , 


(50) 


and the final equation is 


fS(x) - 2 Bf'„(x) + (B 2 + b)Mx) = 0, 


(51) 


where 


B = 


fl(o T 1) 

Now, energy levels are obtained from above and (f|) as 


E k (n) = 9 ^{ 1 


1 Ti Z 

i/7 \ r, — "7 - -- \ b \ —— ~\~ (k~\~Ci)Tl\ 

2 (k + a) 2 (k + n + a ) 2 l 2 


( 52 ) 


The superpotentials are obtained using (|l^), (p]) 
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a k (x) — (k + a)F{x ) + - - 

k + a 


nx) = Bi 


(53) 


where fo(x) is easily obtained in (^Tj). We list the resulting superpotentials in the following 
with a constant d (to be chosen suitably if necessary), and putting 

V = \f\b\ ■ 


(1) b > 0 ; fo(x) = sinrj(x + d) exp (Bx) 


a k (x) 


— (k + a)rj cot rjx + 


9 

k + a 


(54) 


(2) b < 0 ; fo(x) = ( ci sinhrjx + C 2 cosh rjx) exp ( Bx) ; with appropriate choices of 
Ci, c-i, there result two well-known potentials 

a. Eckart potential 

a 



atk(x) = —(k + a)ricothr}x + 


k + 


b. Rosen-Morse potential 


a k (x) = —(k + a)rjtanhr]x + 


k + a 


(56) 


(3) 6 = 0; Coulomb potential ; / 0 (x) = (x + d) exp (Bx) 


. . k + a g 

a k (x) = -+ 


x 


k + a 


(57) 


It is satisfying to see that a single ansatz reflecting the common structure of superpoten¬ 
tials as in (|53|) produces a class of potentials in addition to the starting one, the Coulomb 
potential in this case. The above potentials and the results of our Erst ansatz given in (|23|)- 
(p9|), constitute the whole shape-invariant potentials given in refs. |4|], ||. (There are some 
discrepancies compared with the tables presented in these references.) 

Finally, we want to mention another way of getting superpotentials as following. Since 
it is generally expected that the ground state wavefunctions are products of two factors, or 
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equivalently a k (x) is a sum of two parts, one may try to find new solutions to the difference 
differential equation (jTj) by adding an appropriately chosen second part to an old one which 
is assumed to satisfy it. As a simple ansatz to implement this idea, we take 

Oik{x) = a k F(x) + b k G(x) . (58) 

Here, simple forms are taken for the two parts with constants a k , b k and ^-independent 
functions F(x),G(x) , and a k F(x ) is assumed to solve 0 . Then inserting (|58|) into ( 0 ) 
results in 

(°fe + afc+i) F'(x) + (flfc - ajfc+i) F*(x) = 2 4 a) , ( 59 ) 


(bk + h+i) G'(x) + 2 ( a k b k - a k+1 b k+1 ) F(x)G(x) + (b% - b 2 k+l ) G 2 (x) = 2, (60) 


together with a condition 


(a) . ( b ) 

s k = s y k ’ + sy 


(61) 


Recall that s k should be positive definite, while s ( ' k > , need not be so. In fact, all the 

results obtained so far can also be obtained within this simple ansatz as may be clear from 
the forms of the obtained superpotentials, (2l|), (|38|) , (|53l). We illustrate the procedure by 
obtaining the Coulomb potential : A solution to (|5^) corresponding to s'^ = 0 is directly 
obtained to be proportional to 1/x, a k being determined correspondingly. Next, note that 
([kl) is consistently solved by taking a k b k to be independent of k and G(x) to be a constant, 
say g . Then one may get 


a k (x) = 


k + a 
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x k + a 
One can find h k (x) by taking the difference of 

fo(z) 


V (r) = -l + ( fc + a )( fc + a ~ 1 ) _l 
fcU x + 2x 2 + 2(k + a) 2 


Aa k (x) = 


fo(x) 


- h k (x) , 


(62) 


which can be used to get the fo(x) and h k (x) for the known superpotentials by recalling 
ho(x) = 0. 
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In summary, we have investigated the recurrence relation of Riccati-type differential 
equations in supersymmetric quantum mechanics. Recasting the problem into ([L7]) and 
taking some simple ansdtze on hk(x ), we reproduced the known shape-invariant potentials 
in refs. [Q], 0, and obtained another classes of potentials which look similar but are different 
from the known ones, as illustrated in P2|)-(fT5| ) . It should be stressed that different types 
of potentials can be obtained altogether within a common ansatz on the superpotential as 
shown in the text. By using more intricate ansdtze or introducing some new ideas into the 
formalism, we expect that the method presented in this paper will be helpful for finding out 
more solvable potentials, as well as reproducing the known ones, e.g., Natanzon potentials( 

0, 0, i). 
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